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Abstract. We derive the interaction Hamiltonian of a Laguerre-Gaussian beam with 
a simple atomic system, under the assumption of a small spread of the center of mass 
wave function in comparison with the waist of the Laguerre-Gaussian beam. The center 
, of mass motion of the atomic system is taken into account. Using the properties of 

regular spherical harmonics the internal and center of mass coordinates are separated 
without making any multipolar expansion. Then the selection rules of the internal 
and of the center of mass motion transitions follow immediately. The influence of the 
' winding number of the Laguerre-Gaussian beams on the selection rules and transition 

. probability of the center of mass motion is discussed. 

^ ! PACS numbers: 32.80.Lg, 42.50.Vk 

D^' 1. Introduction 

> ' 

^ I The exploration of the orbital angular momentum (0AM) of the electromagnetic field 
^ I has been impelled by the work of Allen et al. ^\ who pointed out that the Laguerre- 
Gaussian (LG) modes carries a well defined amount of 0AM per unit energy. Since then 
the 0AM properties of the radiation have gained much interest, the concepts of 0AM 
being applied in several fields of optics. Ref. |2I contains a collection of articles devoted 
or in connection to the 0AM of the radiation field. We can take advantage of the 
0AM of the electromagnetic field to study the entanglement in a n-dimensional Hilbert 
space by spontaneous parametric down-conversion [Hj. The interaction of LG modes 
with atoms as point particles has been studied j3] and the transfer of the light 0AM to 
atoms was observed in a four- wave mixing process jS]. Further, by taking into account 
the center of mass (CM) motion of the atomic system, it was shown that radiation field 
endowing 0AM may entangle the internal and external 0AM of the atom [Hj, and the 
field 0AM is imparted between the internal and external motion of the atomic system 
[71 IH]. The aim of this work is to investigate how the selection rules and the probability 
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of transitions of the atomic systems change under the influence of an LG beam, taking 
into account th CM motion of the atomic system. 

2. The atom-radiation system 

We consider here the simplest atomic system, formed by a positive charge +e of mass 
m„ and a negative charge — e of mass nie, both of them spinless. The Hamiltonian of 
the system formed by the atom and the electromagnetic field is given hj H = iJo + -f^int, 
where Ho is the unperturbed Hamiltonian of the atom and i/int is the interaction 
Hamiltonian of the atom with the electromagnetic field, specified ahead in the text. 

The CM coordinate of the atomic system will be given by R = (mgro + mnrn)/mt, 
TTit = me+mn being the total mass, and their relative (internal) coordinate by r = re— Fn, 
where re and rn represent the particles coordinates. Further, we assume that the atom 
state can be written as product of the CM wave function and the internal, i.e. electronic, 
wavefunction T(R, r) = ^j:j(R)\l'(r). The atomic system is considered trapped in a two 
dimensional harmonic potential, and its normalized wavefunction written in cylindrical 
coordinates has the form: 

^r{R) = ^Gnm (—) exp [i{KR, + M<l>)] (1) 

where Gm,m{x) = A/'xl*^lL[^f '(x^) exp(— a;^/2), Lp{x) being the generalized Laguerre 

polynomial, and wr describes the spread of the CM wavefunction. Af = \J2 n~\/n'^\/wR 
is the normalization constant, with = {N ± \M\)/2. The quantum number N gives 
the energy of the harmonic oscillator Er = Tl^ {N + 1) / (wj^mt) , and M gives the amount 
of angular momentum carried by the CM motion Mh. The internal wave function is 
given in terms of spherical harmonics yj"^: 

M/(r) = F„/r)rr(^^,^), (2) 
where Fn^i{r) gives the radial dependence. 

We consider an electromagnetic field, in the paraxial approximation, described by 
a LG beam of waist wq, having no off-axis radial nodes, propagating along the z axis, 
and normalized according to / d^r |i?(r)p = A. The normalization constant A is chosen 
in such way that it is independent of the winding number /. Under the assumption that 
the spread of the CM wave function is smaller in comparison with the waist of the LG 
beam, the electric field reads as [H]: 

En /r I \ I'l 

E(r,t) = ^^ — exp[i{lip + kz - ut)]. (3) 



Next, we write the electromagnetic field in terms of regular solid spherical harmonics 

nr. 

E(r,t) = (-l)^2l'lv^Eo7e(;|(^) exp[i{kz - cut)], (4) 

Wo 

where Eq is the polarization vector of the electric field, and 7^["(r_|_) = C^^r^Y^id = 
7r/2, If) is the regular solid spherical harmonic [HI, with C[" = [Att/ {21 + !)/(/ — m)!/(/ + 
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m)!]^/^. For simplicity, the symmetry axis of the CM motion is considered to be the 
same as the axis of the electromagnetic LG mode. In the subsequent calculations the 
electromagnetic field is treated classically. 

3. The interaction Hamiltonian 

We now focus on the interaction Hamiltonian which plays the central role in revealing 
the influence of the winding number / on system transitions. The structure of the 
electromagnetic field (jH), i.e., in terms of spherical harmonics, is straightway imprinted 
to the interaction Hamiltonian using the Power-Zienau-Wooley scheme |10j : 

H,^t = - / V{r) ■ E(r, t) + h.c, (5) 

where V{r) represent the polarization density. Writing the polarization density as a 
closed integral 

Vir) = ^ e Jr„ - R) dA 5[(r - R - A(r„ - R))] (6) 

a=n,e 

and introducing it in relation (jSJ, the interaction Hamiltonian is expressed as follows 

m 

Hi^t = —v-C dA (mnE(R + A^r, t) + meE(R - A^r, t)] + h.c. = h[^, + (7) 
mt JO I nit ) 

where and refers to the terms given by mnE(R + A^r, t) + h.c. and 

meE(R + A^r, t) + h.c, respectively. 

Using the expression of electromagnetic field in terms of regular spherical 
harmonics, we can separate the CM and electronic coordinates, R and r, respectively, 
by the translation property of regular spherical harmonics jH]: 

^r(x±y) = E i: (±f7^r(x)7^^_];™'(y), (8) 

r=0 m'=-V 

and the electromagnetic field will read as: 




From the properties of the spherical harmonics functions Yj^{9, ip) we know that in the 
case 9 = 7r/2, the functions takes on non-zero values only if / + m is an even number. 
The indices of the spherical harmonics functions corresponding to the internal and CM 
coordinates, in the r.h.s of relation (jH)), are (/', m') and (|/| — l — m'), respectively. This 
implies that when I'+m', also |Z| —/' + /— m' is an even number, so the spherical functions 
takes on non-zero values at the same time and the double sum does not become zero. 
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Introducing the explicit values of the functions Yf^{9 = 7i/2, 0) into equation Q we get: 

E(R + A^r) = 2l'l;MEoexp[iMi?. + A^^)] V V f A^^VV^V""' 
nit ^ m, ;fr"o \ rutWo J Kwq J 

exp(im(^)exp[i(/ - m')^ 

(/' - |m'|)!! (/' + |m'|)!! (|/| - I' - \l - m'\)\\ (|/| - /' + |/ - m'\)\\ ''+"^''^^611 I ) 
where n!! stands for the double factorial. Let us analyze the influence of the double 
factorial (|/| — /' + | / — m'|)!! on the double sum appearing in the above relation. Taking 
into account that \m'\ <l' and V + m' must be even, there are two possible situations: 

(i) / > leads to |/| — /' + |/ — m'| = — /' + m' which takes on either zero or negative 
values. The properties of the double factorial tell us that n\\ = ±oo unless n > 0. 
The case — /' + m' < is ruled out, and we are left with m' = I'. 

(ii) / < leads to |/| — /' + \l — m'\ = —I' — m' also taking on either zero or negative 
values. Performing a similar inference to previous case, the only valid possibility is 
m' = -v. 

Summarizing the conclusion of (jH) and (jHl), the double sum is reduced to a simple sum 
where m' = sgn(Z) I'. The equation (fTUI) becomes 

E(R + A^r) = 2l'l M Eo exp[iA;(i?, + A^^)] V ( A^^V' ( ^) 
JUt ^ m, (r^^ \ rutWo J \wo/ 

exp[isgn(/)/V]exp[isgn(/)(|/| - l')(l>] 
(2/')!! [2(|/| -/')]!! 

In the above relation, one may change the internal coordinate r± to r sin 6 in order 
to give the interaction Hamiltonian in terms of spherical harmonics. Replacing the 
expression of the plane wave by its expansion into spherical harmonics exp(ifcr cos 6') = 



X]^o i^V + ^)jp{kr)Yp{9), where ji,{x) is the spherical Bessel function of the first 

kind of order u, the integration over A is performed and the interaction Hamiltonian 



H^^f^ reads as: 



(2p+l) 

r(2/' + 2) 



oo \l\ 
p=0 Z'=0 o-=0,±l 



I jp+/'(l+sgn{0) 



(i + p + /')r(p+|)r(|/| -/' + !) 



V 2T7i{ / 

( 2±^;p + I, 2±p:^; -{kr^/) + h-c, (12) 



where r ■ Eq was replaced by ryATT/3J2a=o,±i^aYi{6,(p), a being associated with the 
polarization of the electric field, e±i = ±{Ex ± iEy)/\^ and eo = E^. However, in the 
frame of paraxial approximation, the Ez component of the electromagnetic field has a 
small contribution to the total electric field. The latter one can be given, to a good 
approximation, in terms of the transverse electric field 
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In relation (jl2j) . since p runs from zero to infinity and the hypergeometric function 
1F2 consists of an infinitely series, it follows that all multipole transitions are excited. 
Nevertheless, the factor r^'+P+'^ allows us to consider the first and the second transitions, 
while keeping only the first term, i.e., the constant one, of the hypergeometric function 
1F2: I' + p + 1 = 1 gives the dipole transition, I' + p + 1 = 2 gives the quadrupole 
transition. We note that in the dipole approximation, there is no transfer of the field 
0AM to the electronic motion, our result being consistent with the concluding remark 
ofRef. m. 

4. Selection rules and transition probabilities 

The selection rules of a certain transition are obtained by the transition matrix element 
between the final and initial state A^i^/ = {'^ f\Hinf^\Ti). Putting together the 
expression of the interaction Hamiltonian (fT^ and of the atomic wave function (|T|- 
(|21), the transition matrix element is written as: 



In relation (fT^ the Dirac delta function and Kronecker delta symbol express the 
conservation of the axial component of the CM momentum and the conservation of 
CM 0AM, respectively. The bracket symbols designate the integral over the radial 
components R± and r, and over the internal solid angle element d{cos6)dip, respectively. 
The transition probabilities are proportional to the integrals over the radial coordinates 
R± and r, being directly influenced by the absolute value of the winding number I. 

The selection rules which governs the electronic-type transitions are given by the 
integral involving the spherical harmonics in relation (fT!^ . Using the addition theorem 
of angular momentum, the spherical harmonics can be coupled and the final result 
expressed in terms of Clebsch-Gordan coefficients. The parity selection rules can 
be easily determined taking into account the properties of spherical harmonics, the 
transition being allowed only when Ij + li + l'+p + l is an even number. From 
Table [T] one can see that the selection rules of the dipole and quadrupole transition 
are unchanged, as expected. The z component of the angular momentum of the entire 
system is conserved, i.e., the relation Am + AM = / + a is always fulfilled pi. Here, 
since we are in the paraxial approximation, the polarization of the electromagnetic field 
determines changes only in the internal magnetic quantum number m of the atom, but 




" [r(2Z' + 2)J (/'+p + i)r(p + 3/2)r(|/| -/' + !) 



(13) 
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Transitions driven by field 0AM 
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Table 1. Selection rules of transitions induced by an LG beam: (A^, Am) and AM are 
the changes of the electronic-type quantum numbers and of the CM 0AM, respectively. 
The numbers under column of Am correspond to different polarizations a of the 
electromagnetic field. 



beyond paraxial approximation this is no longer valid ^2]- In relation as /' runs 
from zero to |/| the 0AM of the electromagnetic field is imparted to the CM motion 
and to the electronic-type motion. The maximum value of the transferred angular 
momentum either to the CM motion or to the electronic-type motion does not exceed 
|/|. This feature is a pecuharity due to the radial dependence of the electromagnetic 
field, assumed in relation Q; however, in general, no bound on the value of the angular 
momentum transferred is set |7| . Although the selection rules are not directly governed 
by the winding number / of the LG beam, they are sensitive to the sign of / as one can 
see from Tabled in the case of the quadrupole transitions with p = and I' = 1. 

The influence of the winding number / on the transition probability can be 
analyzed by performing the radial integral over the CM coordinate R± in formula (|T^. 
Considering the initial and final CM states of the form (^, the integral is given by ^ 
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with Mf = Mi + sgn 
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nl 



n 



f 



1/2 



+ 1: 1 + \Mi 



\Mi\~Nf + \l\-l' 



+ i;i), (14) 



-n 



, Depending on the sign of the product / Mi, the relation 
(fTl]) may be further simplified. 

Next, we will focus on the CM transitions due to the photon absorption in the 
dipole interaction, p = I' = 0, and in the quadrupole interaction, with one 0AM unit 
transferred to the electronic motion, p = and /' = 1. The dependence of the CM 
transition probability Pcm on /, /' and Wr/wq is the form 

2 



Pcm oc 



/'(l+sgn(0) 



/ 4 



;i + /')r(|/| -/' + !) \kwo 



r(2/' + 2) 



{GNf,Mf I 



kRj 



\i\-i' 



\Gn„au) 

(15) 
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The selection rules of CM transitions are given by (i) the properties of the ID- 
harmonic oscillator (through the radial integral (fT^ ) and (ii) angular momentum 
conservation law (through the variable $), which fixes the quantum numbers Nf and 
Mf. In figure HI we plot the transition probability Pcm as function of Nf for / = 2 (left) 
and I = 3 (right). In figure^, the CM selection rules are AA^ = Nf — Ni = 0, ±2 (with 
the electronic transition of the dipole type, /' = 0), and AA^ = ±1 (with the electronic 
transition of the quadrupole type, /' = 1). In figure^, for / = 3, the number of the 
final CM states which can be addressed is larger in comparison with figure^, due to 
the multipole CM transition |/| — /' = 2, 3. 



0.01 - 
le-04 - 

le-06 - n 
„o le-08 - 
le-10 - 
le-12 - 

le-14 ' ' U ^ ' 

2 3 4 5 6 

Nf 

(a) (b) 

Figure 1. The CM transition probability Pcm as function of the energy quantum 
number Nf, for N, = 6, Mi = 0, wr/wq = 10"'', I = 2 (left) and / = 3 (right). 
One can see how the interaction type, i.e. dipole or quadrupole, changes with I for 
the final states. The white and black boxes correspond to CM transitions due to 
the electronic dipole {I' = 0) and quadrupole (/' = 1) interaction, respectively. The 
transition probabilities are normalized to the the value corresponding to Nf = 6 in the 
left figure. 

The quantum numbers of energy and of the angular momentum M, appearing in 
relation must be of the same parity (this comes from the fact that the wavefunction 
describes two uncoupled ID-harmonic oscillators with the quantum numbers and 
Uy, and from here follows N = + Uy and M = — Uy). When one considers 
only the interaction in the leading orders, the difference in the quantum number Mf 
will be — Mj^^ = ±1, where (d) and (q) designates the dipole and quadrupole 
interaction, respectively. Therefore the energy number of final states reached by dipole 
interaction will have opposite parity in comparison with that one of final states reached 
by quadrupole interaction, see also figure ^ It follows that we can not address using 
only one LG beam a specific state by both dipole and quadrupole interaction. To enable 
both interaction types one may use superpositions of LG beams having different winding 
numbers /, e.g. by suitable placed non-axial vortices [Hj. 

From figure Q one sees that the CM transitions occurring due to the electronic 
dipole interaction are about 10^ times weaker than those occurring due to the electronic 
quadrupole interaction. However, when one calculates the overall transition probability, 
i.e., electronic and CM, the difference between these two transition types will decrease 
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to about 10^ times, because the electronic quadrupole transitions are about 10® times 
more unlikely that the electronic dipole transitions. 

In figure El we plot the CM transition probability as a function of the winding 
number / and of the ratio wh/wq, between the initial Ni = 6, Mj = and final Nf = 12 
state. The transitions are driven by the dipole (quadrupole) interaction for / even 
(odd). Introducing the expression of the matrix element ()14|) into relation (jlSj) the 




LG winding number 1 



Figure 2. The CM transition probability Pcm versus the winding number I and the 
ratio wh/wq. The initial state is characterized by TVi = 6 and Mi = 0, and the final 
state by Nf = 12 and Mf — sgn(/)(|^| — I'). The values were normalized to Pcm for 
I = 0. 

dependence of the transition probability Pcm on the ratio wr/wq will be of the form 
Pcm oc (wij/wo)'''^''- This behavior characterizes the approximation wji much smaller 
than Wo, and one may compare it with the dependence obtained from formulas derived 
in Ref. |7] in order to establish the validity range of the approximation. 

5. Conclusions 

The selection rules and the matrix elements ()13|) shows that dipole transitions affect only 
the CM transition. On the contrary, for quadrupole transitions there is an exchange 
of one unit of angular momentum between the LG beam and the electronic states, 
and no limitation exists in the exchange between the CM and angular momentum 
of the beam [see first bracket in (fT!?|) ]. In the case of an atomic system trapped 
in an harmonic potential the dependence of the CM transition (selection rules and 
transition probability) on the LG winding number I is discussed. There are systems 
[T3] like silver nanoparticle where the quadrupole transition are even more intense than 
dipolar transitions, so that the difference of about two orders of magnitude between the 
electronic dipole and quadrupole transition can be decreased. In this case we suggest 
to perform absorption experiment using LG beam to test the dependence of quadrupole 
transition on /. Very likely, due to the correspondence of quadrupole transformation 
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rules and polarizability tensor, the dependence of CM motion on /, found in this paper 
and in Ref. [7j, would affect Raman and Brillouin transition intensity and the relative 
selection rules. This would further increase the interest in using LG beam for novel 
spectroscopy studies. The application of LG beams to Raman and Brillouin transitions 
will be investigated in a further work. We conclude pointing out that even if the internal 
transition in quadrupole approximation is limited to one unit of 0AM, the dependence of 
CM transition on / is of interest in non resonant transition, where the angular momentum 
of the LG beam, can be used to modulate the transition probability. 
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